Orchestrating an NMR quantum computation: the N = 3 Deutsch-Jozsa algorithm. 
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A detailed description of the development of a three qubit NMR realization of the Deutsch- 
Jozsa algorithm [Collins et.al., Phys. Rev. A 62 022304 (2000)] is provided. The theoretical and 
experimental techniques used for the reduction of the algorithm's evolution steps into a sequence 
of NMR pulses are discussed at length. This includes the description of general pulse sequence 
compilation techniques, various schemes for indirectly coupled gate realizations, experimental pulse 
| parameterization techniques and bookkeeping methods for pulse phases. 

(N 

^\ I. INTRODUCTION 

, 

Nuclear magnetic resonance (NMR) spectroscopy has become well established as a tool for exploratory experimental 
investigations of quantum computation. Fundamental gates, algorithms, error correction and other issues in quantum 
information theory have been demonstrated at the level of a few qubits using room temperature, solution state NMR 
systems, [fl]— 126| 

Recently we successfully implemented the N = 3 case of the Deutsch-Jozsa algorithm using room temperature 
solution state NMR. |l2| This implementation was characterized by the choice of qubits, namely the 13 C nuclei of 
alanine, and the resulting methods for processing the quantum information. The homonuclear nature of the qubits 
(henceforth qubit and spin are used interchangeably) demands qubit selective pulses throughout the realization, thus 
(/"■) ' greatly complicating the experimental effort in comparison to a heteronuclear system. Furthermore the arrangement 
of couplings between qubits in alanine is such that the only reliable method for processing quantum information is 
along a chain of qubits, whose pattern of couplings is effectively A — B — C (i.e. the A — C coupling is too weak). This 
demands "indirect" gate realizations using swapping techniques via intermediate spins p7{ and is more difficult than 
■ "direct" realizations, possible in the fully coupled case. Although these restrictions are artificial at the level of a few 
i[ qubits there is reason to believe that they will become important as the number of qubits increases p6[ ]. Therefore 
^Jj i the ability to implement the associated information processing techniques will become crucial; hence our choice of 
this restrictive qubit system. 

Even at the low level of three qubits the resulting experiments were sufficiently complicated to motivate the 
development and use of a host of simplification techniques and experimental methods to construct the appropriate 
pulse sequence. Much of the literature bypasses or only briefly describes such details, which will become increasingly 
important as the number of qubits in NMR quantum computation implementations grows. The purpose of this 
article is to describe in detail the construction of an NMR realization from the theoretical algorithm through to 
the final spectrometer output and the algorithm result. Section || briefly describes the Deutsch-Jozsa problem and 



* 

algorithm and the initial decomposition of the algorithm's evolution stage into fundamental quantum gates. Section III 
covers the theoretical development of initialization schemes, pulse sequences and readout schemes. General pulse 
sequence compilation techniques and methods for constructing indirectly coupled gate realizations are presented here. 



Section [V describes the experimental details and includes pulse parameterization, bookkeeping methods for pulse 



phases, construction of appropriate evolution delays, results and possible sources of error. 



II. DEUTSCH-JOZSA ALGORITHM 

The Deutsch problem [ p8|j29| addresses a global property of certain binary valued functions / : {0, 1} N — > {0, 1}. 
In particular, a function is balanced if it returns as many times as 1 after evaluation over its entire range. Given any 
function which is either balanced or constant the problem is to determine its type. Classical algorithms which rely on 
repeatedly evaluating the function for various arguments require a number of evaluations which grows exponentially 
with N to answer the problem with certainty J2S| ] . However, a quantum algorithm requires a single evaluation off 
and solves the problem with certainty |2g| |30(| . In its simplest rendition the algorithm needs an N qubit register |3(]] 
and follows the scheme illustrated in Fig. |U 
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FIG. 1. Deutsch-Jozsa algorithm. Initialization to the state: \x) = \xn-i) ■■■ \%o) where Xi € {0,1}. J? n (90) rotates each 



qubit through 90° about n = cosc/ix + sin^y. For a single qubit R n (9) 
/-controlled gate, Uf \x) — < -^ !{x) 
conclusively. 
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where n.cr ~ n x a x + n y a y + n z a z . The 
( — \x), evaluates /. The expectation value of \x) {x\ on the output state answers the problem 



The ease of implementation relative to other quantum algorithms has made the Deutsch-Jozsa algorithm a natural 
candidate for experimental NMR demonstrations of features of quantum computation. The N — 1 case, using the 
Cleve version [ p9[ , was the first quantum algorithm implemented experimentally |l],|| . This was followed by the N = 2 



case 



11 



However, the algorithm is uninteresting from a quantum computational perspective for N < 3, in which 
case the essential operations could be implemented classically |}0|. Recently the modified version (see Fig. 0) of the 
N = 3 algorithm has been implemented for various balanced functions |I^-|l4[]. Finally, limited balanced function 
cases at the N = 4 level have been tested using the Cleve version |l(]] . 



A. Gate Decomposition 

The heart of the Deutsch-Jozsa algorithm is the function evaluation step, 

U f \x):=(-lf x) \x), (1) 

and the bulk of the effort in any experiment will be devoted to its realization. This is initiated by decomposing Uf 
into a product of fundamental one and two qubit operations which are amenable to physical realization. In general 
single bit rotations (one qubit operations) and controlled-NOT gates (two qubit operations) are satisfactory. |HJ In 
this scheme, the hope is to be able to devise physical realizations for the fundamental gates and concatenate these 
to recreate U f . The definition of Uf, Eq. (|l|) , does not immediately provide such a decomposition. One approach 
for arriving at a fundamental gate decomposition begins by representing any admissible function via power series 
expansions in the argument bits, Xi where i = 0, 1, 2. Such expansions terminate since xf — Xi for Xi G {0, 1}, leaving 
at most a cubic term for N = 3. However, it can be shown that for balanced and constant functions the cubic term 
disappears. Thus for N = 3 any admissible function can be expanded as: 

2 2 

f(x 2 ,xi,x ) = (J) a ij x l x j ^a i x i ® a (2) 

i>j>0 i=0 

where addition is modulo 2 and , a^, a € {0, 1}. Therefore 

2 2 

Uf\x2x 1X0 )= n (-lr^n^^'i-^iwo}. (3) 

i>j>0 k=0 

The constant term merely provides an irrelevant overall phase factor and can be ignored. The quadratic and linear 
terms motivate the definitions of the quadratic term gates and linear term gates: 

U% VAD \x) :=(-ir^» and (4) 
UL\x) :=(-lp\x), (5) 

leaving 

2 2 

Uf= n (^Ao) a,j n(^™) Qfc (6) 

i>j>0 fe=0 
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where the order of operation reads from right to left. However, quadratic and linear term gates all commute and any 
rearrangement is permitted in this decomposition. Decompositions of these in terms of fundamental gates are 

Ut m = i?|(180°) and (7) 

where J/cnot is a controllcd-NOT gate with control i and target j and R^O) is the operator representation of a 
rotation of spin i about axis n through angle 9. 

B. Function Classification 

There are f njv— i) balanced functions with N bit arguments and for ./V — 3 this gives 70 distinct possibilities. However, 

certain distinct balanced functions are such that the operations required to realize Uf are related by permutations of 
the argument qubits. For example, permuting the qubit indices using the rule 1 <-» 2 transforms the operation required 
to realize f{x2,xi,x<P) — X2 to that for g(x^,xi,xo) — x\. The physical difference between the realization for these 
cases is minor in comparison to that for a case which is unrelated in this fashion; an example is h{x% , X\, Xq) — x-i © x\ . 
Thus it makes sense to classify the balanced functions according to the similarity of their power series expansions 
under permutations of argument indices; this is mirrored by an equivalent classification of the operations required for 
function evaluation under qubit permutations. Within each class the arrangement of linear and quadratic term gates 
differ only by permutations of their qubit indices and this reduces the number of essentially different experimental 
realizations needed to demonstrate the algorithm. It emerges that there are ten classes of balanced functions; a 
representative of each is provided in Table |. All possible quadratic term gates are required for the classes represented 
by fg and fio and all possible linear term gates are required for ^3. Thus realization of the algorithm for one 
representative from each class requires the construction of all possible linear and quadratic term gates, which could 
be used for realization of any other case. The experiment will be limited to one function from each balanced class 
plus the constant functions. 
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TABLE I. Representatives of each class of admissible functions (N = 3) and the corresponding density operator after the 
function evaluation step. 



III. PULSE SEQUENCE DEVELOPMENT: THEORY 

NMR spectroscopy of spin i nuclei of appropriate molecules in solution offers a readily accessible experimental 
approach to quantum computing. Any molecule containing three distinguishable, coupled spin ^ nuclei in an 

external magnetic field provides the three qubits needed to solve the N — 3 Deutsch problem. To a good approximation 
the Hamiltonian for a room temperature, solution-state sample is 
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2 2 

#=EiH+f E ( 9 ) 

8=0 i>j>0 

where are the Zeeman frequencies, Jy the scalar coupling constants, and a\ Pauli operators f3^] . Superscripts 
label the spins and identify them with the corresponding argument bits. 



A. Gate construction and compilation 

Equations (^J) - (|J) provide a method for constructing a real pulse sequence for Uf in terms of an appropriate 
sequence of single bit rotations and controlled-NOT gates. Pulse sequences for the latter then enable a practical 
realization of the algorithm. The resulting sequences can be displayed, following conventional quantum mechanics 
formalism, as a sequence of operators with order of application read from right to left, U m • ■ ■ U% U\. An alternative, 

more familiar to NMR users, is a sequence of pulses, \U{\ — [U2] •— [U m ], with order of application read from left to 

right. In the room temperature, solution-state NMR paradigm two types of operations are necessary and sufficient for 
gate construction. These are spin selective rotations caused by appropriate externally applied RF pulses, and delays 
during which the system evolves under only one of the Hamiltonian coupling terms. The process of pulse sequence 
construction is initially simplified by considering spin selective rotations which are idealized in the sense that they are 
instantaneous and yet perfectly selective. An ideal rotation of spin i about axis n through angle 9 will be represented 
by the pulse sequence term [9f n , corresponding to the operator R\(0). A delay or period of free evolution of duration 
t under the scalar coupling term ^ Jija\a\ (no sum over the indices) alone will be represented by 

UtL,(t) ■= e-'W'* (10) 

or the pulse sequence term [t] lJ . This operation is unnatural in the sense that some of the scalar coupling terms are 
absent. In practice it is more convenient to build this from periods of free evolution of duration t under all scalar 
couplings, one of which will be denoted 

UlZUt)-=e~ iT '^^ JskHait (11) 

or [t] tot . 

Single qubit gates, equivalent to qubit selective rotations, can be implemented by applying spin selective RF pulses 
of the appropriate shape, power and duration J32|. Controlled-NOT gates can be realized using the rotation and delay 
method |],|||. There are two idealized decompositions (i.e. in terms of idealized rotations): 

tf&or = #.(90°) R^(90°) i4*(90°) ^ CAL (l/2J y ) J^(90°) (12) 
or in pulse sequence terms: 

[Ug OT ] = [90T ±a - [U^Jijf - [90°]L ~ [ 90 °]L - [90°]i 2 (13) 

where the axes of rotation are given by either all the upper signs or else all the lower signs. 

In order to implement any algorithm, it suffices to simply concatenate constituent pulse sequences for the required 
fundamental gates. However, the difficulty of experimental implementation may be reduced by various simplification 
techniques, which can be implemented efficiently at the algorithm design level via rudimentary, general schemes, 
applicable in any physical realization. These use 

RU(0) Rktt) RUO) = Rl(J>) (14) 

where 

f = mcos# — n x msinS + (1 — cos6>)(n • m)n (15) 
and which is easily verified by expansion and multiplication. This provides two avenues for simplification: 

1. The pulse sequence can be rearranged so that all rotations about z follow all other operations. This is possible 
in the context of the linear and quadratic term gates because such rotations commute with periods of evolution 
under H and because they can be interchanged with rotations about axes in the transverse plane by applying 
Eq. (here illustrated for qubit i, but applicable for any i): 
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#cosax+smay(0) R li d ) = R W) Kos (o!-0)x+sin (o!-0)y(# ( 16 ) 

The result is that, all rotations about z on qubit i can be replaced, after modification of the succeeding rotations 
on qubit i, by a single rotation about z on qubit i immediately prior to the measurement step. In the context 
of solution state NMR spectroscopy this can be implemented via a qubit dependent phase shift applied to the 
acquired data. The need to actively implement i rotations is thereby eliminated. 

2. Successive 90° rotations (on one qubit and with no intervening period of evolution under couplings involving 
this qubit) about orthogonal axes in the xy plane, rh and n. can be replaced by a single rotation about an axis 
in the xy plane followed by a rotation about 2. Again Eq. ( |14| ) prescribes: 

£1(90°) J&(90°) = Ri xA (90°) £1(90°) (17) 

and n x m is oriented along z. The z rotation can be interchanged with all succeeding operations according to 
the previous scheme. This reduces the number of rotations requiring active implementation. 

The number of steps required to carry out these compilation procedures scales polynomially with the number of 
gates and thus does not adversely affect the computational complexity of the algorithm. Applying these rules to the 
concatenated sequences for the Deutsch-Jozsa algorithm yields: 

1. The linear term gates and other single qubit rotations about z can be interchanged with the remaining gates 
(with appropriate modifications to their axes of rotation) and absorbed in the final qubit-dependent phase shifts. 
These can be implemented passively via mathematical transformations on the acquired data and need not be 
included in the actual physical information processing stage. The fact that this gives a single manipulation per 
qubit means that this strategy does not change the computational complexity of the problem. 

2. The remaining gates in the sequence are either single qubit rotations through 90° about the x or y axes and 
scalar coupling gates. For any qubit, the single qubit rotations which appear between successive scalar coupling 
gates involving this qubit can be reduced to at most one rotation about an axis in the xy plane and at most one 
rotation about z, which can be interchanged and absorbed as a phase shift. 

The result of compilation is a sequence of scalar coupling gates interspersed with single qubit rotations about axes in 
the xy plane and terminated with one qubit dependent rotation about z per qubit. Thus the sequence of operations 
involving qubit i (ignoring operations which only involve the remaining qubits) has the form, 

WW, a \0\V) = WM) Rl osa , - +sinQl J9 m+ i) tfe(4) • • • 

m+l tjt+1* 7 

' ' ' ^rasa'x+sma'y^) ^scalW) Kob a{ X+sin a{f (#0 ( 18 ) 

where m is the number of scalar coupling gates in the sequence, and a 1 := {a^+ii • • • i a \}i ^ := {^m+i> • • • > ^1} an< ^ 
V := {t l m , . . . , t\} are arrays of parameters which are determined by the algorithm and pulse sequence compilation. 

The distillation, via the compilation scheme, of all qubit selective rotations about i to one rotation per qubit 
immediately prior to measurement and the possibility of implementation via post-acquisition mathematical operations 
motivates the division of the pulse sequence into an active section, 

i>: ctiv >\0\t 4 ) := £L< +lX+sin < +iy 0Wi) • ■ • uii\M nL^+^afrVi), (19) 

and a passive section, 

uUJtf) ■■= Rim, (20) 

giving 

u\<p\ a \e\e) = ui ssiv M) eSt'). (21) 

In an NMR realization the active section has to be implemented by externally applied fields and delays for internal 
Hamiltonian evolution. This requires an experiment in which the appropriate pulse parameters (corresponding to a 1 
and 6 l ) and delays (corresponding to V) are carefully adjusted so as to achieve the correct evolution. In contrast 
the passive section is best implemented mathematically on the acquired data using the spectrometer's phase shifting 
capability. 
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An interesting and useful feature of the active section of the pulse sequence is that the absolute orientations of its 
constituent rotations' axes are less important than their relative orientations provided that the initial density operator 
satisfies certain general conditions. This results from Eq. (Q) which gives 

-^cosax+sinayW = ^"-z(V? 1 ) -^cos (a-^ ! )x+sin (a- V ')y (^) R\vP ) (22) 

and consequently 

u\0\ a \8\t>) = ui^M - v*) uL^W - v'^t') W). (23) 

where a 1 — cp l := {<x l m+1 — (p l , . . . , a — tp 1 }. The initial rotation, R' l -(ip l ), can be omitted whenever the initial state 
of the system is a mixture of eigenstates of I z operators. Examples include the thermal equilibrium state for weakly 
coupled nuclei and pure (or the equivalent pseudopure) states such as \xn-i ■ ■ ■ Xo)- Henceforth we shall consider this 
case, so that 

U\<P\c*\e\e) = L>^ aivo (0 l - <?) - V\0\t*) (24) 

and which demonstrates that, for qubit i, only the relative orientations of the axes of the rotations on qubit i are 
important. Thus the axis of the initial rotation for each qubit can be chosen arbitrarily with no need to correlate 
these between the different qubits. Thereafter the axes of the succeeding rotations are fixed. The freedom of choice 
of the axis of the initial rotation is offset by the necessity for a compensatory adjustment in phase of the passive 
section rotation although later it will emerge that this is unimportant for the N = 3 Deutsch-Jozsa algorithm. This 
simplification is akin to working in co-rotating frames (one for each spin) which do not coincide at the point of the 
first spin selective rotation and Eq. (p3) demonstrates that this is legitimate. 

B. Linear and quadratic term gate construction 

For linear term gates the resulting idealized pulse sequence is 

[UL] = [180°]; . (25) 

The function evaluation gate can always be arranged so that the linear term gates appear after the quadratic term 
gates. This adds f80° to the passive section z rotation for spin i and can be implemented as part of the post-acquisition 
phase shift. This motivates a further classification of admissible functions based on the similarity of the active sections 
of the evaluation step pulse sequence (i.e. on the similarity of the quadratic terms in /). Ignoring all linear terms and 
distinguishing only between arrangements of quadratic terms, the following classes emerge: /2) /a}, {fi, /s, fe}, 
{/V,/g} and {/g,/io}- Within a given class the pulse sequences differ only by combinations of 180° post-acquisition 
phase shifts; their active sections are identical. There are essentially four different experiments to perform to implement 
N = 3 Deutsch-Jozsa algorithm. 

Apart from the initial rotations, the only other contributions to the active sections of the pulse sequences come 
from the quadratic term gates. Substitution of a controlled-NOT gate pulse sequence, Eq. (|l2|), into the quadratic 
term gate decomposition, Eq. (||), followed by compilation yields the idealized pulse sequence 

Ug AD = #_,(90°) RL(90°) l^ CAL (l/2J y -) (26) 

or 

[Ug AD ] = [1/2J«]« - [90°Y_ Z - [90°]i* • (27) 

The success of a literal translation of this sequence into a refocused delay followed by relevant phase shifts depends 
on the ratio of \Jy | to the fastest relaxation rate for the spins which represent the qubits. If this ratio is large enough 
then such a directly coupled realization is feasible. 

On the other hand, if \Jij \ is comparable to or less than any of the relaxation rates the directly coupled realization 
will result in unacceptable errors. However, provided that a network of substantially stronger couplings connects i 
and j it may be possible to implement an indirectly coupled realization by swapping information along the coupling 
pathway. Suppose, for example, that each of i and j are coupled to k appropriately strongly. Then the SWAP 
operation U^, defined on the computational basis as 

t^swAP \ X N ■ ■ ■ Xi . . . x k ■ ■ ■ x ) := \x N . . . Xk ■ ■ ■ Xi . . . x ) (28) 

and extended linearly to all states, allows || information to be exchanged via k. A decomposition into controlled-NOT 
gates is given in Fig. |2[ Two possibles strategies for the construction of U^ CAIj (l/2Jij) are apparent: 



G 



FIG. 2. Decomposition of the SWAP operation, E/swapi into controlled-NOT gates. 

1. Use evolution under the coupling between k and j and the swapping technique to give 

^scAiXl/'^j) — U SWAP J7| CAL (l/2 Jjk) U SWAP . (29) 

One possible compiled pulse sequence is provided in Fig. [| Others are available by mixing alternative construc- 
tions for the constituent controlled-NOT gates. 
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FIG. 3. Pulse sequence for indirectly coupled realization of Ul J CAJj (l/2Jy) based on swapping to a scalar coupling gate. Lines 
connecting qubits i and k represent U$% AIj (l/2Jjfc). Squares represent single qubit rotations. 



2. Use an indirectly coupled controlled-NOT gate between j and k to construct a controlled-NOT gate between i 



and j, which in turn gives U S ca^0- /2</y)- Thus 



fpi = lP k JJ ik fP k 

u CNOT U SWAP u CNOT ^ SWAP 



(30) 



and 



^ CAL (l/2J. y ) = # -(90°) [7|£ AP L>f NOT U»„ 4(90°) 



(31) 



and simplification gives the sequence of Fig. |J. 
The relevant quadratic term gate can then be constructed by appending the appropriate rotations about z. 
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FIG. 4. Pulse sequence for indirectly coupled realization of (7scal (1/2 Jij) based on swapping to a controlled-NOT gate. 
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C. Management of Hamiltonian evolution 



The pulse sequences of Eqs. (|T^) and or Figs. ^ and ^ assume that the internal Hamiltonian for the molecule, 
Eq. (||), is only active is during the delays which appear in the sequence. Furthermore during such delays it is 
assumed that only one scalar coupling term is active. These assumptions fail in two ways. First, externally applied 
spin selective RF pulses have a finite duration, during which the Hamiltonian is active. Second, during any delays, 
which are implemented in practice by allowing the spins to evolve freely for an appropriate duration, all the terms of 
the Hamiltonian are active. 

During the externally applied pulses this is addressed by neglecting evolution under the coupling terms, which are 
typically small in comparison to the external field strengths, and noting that evolution under the Zeeman terms can 
be viewed as rotations about the z axis. In this approximation the effect of Hamiltonian evolution has been to include 
unwanted qubit-dependent z rotations. Fortunately, Eq. (|l6|) provides a rule for interchanging any such z rotation 
with succeeding rotations about transverse axes for the same qubit and gives the necessary phase correction for these. 
Additionally they commute with any scalar coupling terms. Thus the Zeeman evolution terms can be interchanged 
so that they effectively appear at the end of the algorithm pulse sequence. This procedure amounts to working in the 
co-rotating frames, which essentially remove the Zeeman terms from H. 

The same technique effectively removes the Zeeman terms during delays for scalar coupling evolution. The remaining 
difficulty is to "turn off" evolution due to other scalar couplings during the delay. This is accomplished via standard 
refocusing schemes p^-|36| , which are based on periods of free evolution under all scalar coupling terms interspersed 
with qubit selective 180° rotations. For example, the sequence 

i4(180°) e ~ l " T/2 i4(180°) e -' lflT/2 (32) 

where rh = m x x + m y y, results in evolution of duration T under all scalar couplings except those involving spin j 
together with the identity operation on spin i (ignoring the Zeeman evolutions as discussed previously). 



D. Initialization and readout 



A typical room temperature solution state NMR sample consists of an ensemble of non-interacting identical 
molecules. Although the nuclear spin state of any individual constituent molecule is generally inaccessible, it is 
possible to control and acquire information about the ensemble's average state, typically described in terms of the 
system's density operator. Fortunately, it has been possible to modify conventional quantum computation schemes ac- 
cordingly. In particular, initialization schemes, relying on the concept of pseudopure states have been developed M ||] 
and tomography schemes allow for the reconstruction of final state density operators [Q. The algorithm stage is 
implemented by subjecting the ensemble to a pulse sequence derived from the concatenation of the pulse sequences 
for the algorithm's constituent fundamental gates. 

Bulk state preparation and tomography schemes typically require efforts, in the form of additional pulses or repeated 
experiments, beyond those demanded by the algorithms. However, it emerges that, for N < 3, it is possible to apply 
the evolution stage (see Fig. |l|) directly to a thermal equilibrium initial state and successfully solve the Deutsch 
problem with an expectation value measurement |3"|]33]|. Analysis of the system state after the function evaluation 
step, which is equivalent to running the algorithm followed by the readout Ra(90°), demonstrates this for N — 3 (the 
reasoning is similar for N < 2). Using the product operator formalism, the deviation part of the thermal equilibrium 



equilibrium density operator for a weakly coupled homonuclear NMR system is proportional to 32 

p tYl = P z +il+P z . (33) 
Setting the initial rotation's axis to n = — y results in 

p th ^ 0O) -il - n - 4° 04) 

^ Pf (35) 
which are listed in Table ffl. The spectrometer measures the components of the magnetization along the x and y axes, 

2 

(M + fo>)) (*) := E Tr {% + d l) e-^pj"^) , (36) 
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where p is the density operator immediately prior to acquisition. Finally the spectrum resulting from a time domain 
Fourier transform of (-M+ (</?)) (t) is calculated and displayed. The acquisition phase, ip, describes any overall phase 
which is independent of the pulse sequence applied to the sample. Typically it depends on factors such as the spec- 
trometer's construction, probe construction or receiver electronics and may be difficult to determine. The appearance 
of any spectrum depends on ip, although this can be adjusted arbitrarily by applying phase shifts to the acquired 
data. In conventional NMR applications the lack of knowledge of ip does not usually impose any impediments. For 
example, post-acquisition phase shifts can be applied to the data so a line of interest appears in absorption mode. 
However, for quantum computation purposes tp could be important because the only possible method of extracting 
information from p, and hence results from the computation, is via the acquired spectrum. Fortunately, the fact that 
p is constant for different pulse sequences means that it is completely irrelevant when comparing two spectra provided 
that the physical setup (except for the pulse sequences) used to acquire these is unaltered. In practice such compar- 
isons are facilitated by adjusting post-acquisition phase shifts. It follows that the algorithm will have to be answered 
by comparing spectra. Signal acquisition immediately after application of Uf and with no additional readout pulses 
provides a spectrum, the f -spectrum, whose form depends on /. This will be compared to a fiducial spectrum which 
is obtained in the same way but with Uf replaced by /; here the system's pre-acquisition state is pad = —1% — I x — 1%, 
identical to the /-spectrum for a constant function. The pertinent features of the resulting spectra can be deduced 
from the density operator prior to acquisition and depend on the presence of scalar coupling terms. The analysis 
is simplest in the case where all scalar coupling strengths are significantly larger than the line widths. Here each 



a) 



b) 
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FIG. 5. Multiplet forms: a) in-phase quartet arising from P x , b) antiphase quartet arising from P z Ii, and c) doubly antiphase 
quartet arising from P Z P X I Z . 

multiplet consists of four distinct lines. The form of the spin j multiplet (i.e. the lines' relative intensities and phases) 
depends on the density operator term containing a factor of V x (see Fig. ||). A single factor of I x yields a multiplet 
with lines of equal phases and intensities. P Z I 3 X yields an multiplet in which the pairs of lines separated by frequency 
Jij have opposite phases but equal intensities (an antiphase multiplet). P Z I X I Z yields a multiplet in which pairs of 
lines separated by a Jij have opposite phases while simultaneously those separated by Jkj have opposite phases and 
equal intensities (a doubly antiphase multiplet). Multiplication of a product operator term by —1 corresponds to a 
phase shift of 180° for the corresponding multiplet. Now consider the case where one or more of the scalar couplings 
is small in comparison with the linewidths. Suppose, for example, that Jij is small relative to the linewidths for spin 
j. It follows that the spin j multiplet's constituent lines which would otherwise be separated by Jij will overlap. The 
resulting lineshapes can be predicted by adding two appropriately overlapping Lorentzian lines. In the case of lines 
of opposite phase separated by Jij partial cancellation will result. In the extreme case, Jij = 0, the cancellation is 
complete and these lines disappear. Thus the multiplet arising from any product operator term containing a factor 
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of P z li disappears if Jjj = 0. 

The quantum logic for this algorithm demands that for each qubit is coupled to at least one other qubit with a 
scalar coupling much larger than any relaxation rates (and hence linewidths) . For a three qubit molecule this implies 
that only one scalar coupling term can be negligible. In this case Table | indicates that at for fj to /io at least one 
antiphase or doubly antiphase multiplet survives. For to fe the multiplet for spin survives but its lines are all 
out of phase with those of the fiducial spectrum. For f\ to f% the same applies to spin 2. Similar results apply for 
other representatives of the function categories but with the qubits permuted appropriately. On the other hand the 
spectrum for the constant functions is identical to the fiducial spectrum. Thus we arrive at an algorithm answer 
criterion for applicable functions: 

(i) / is constant if and only if the f -spectrum is identical to the fiducial spectrum and 

(ii) / is balanced if and only if there is a n phase difference between at least one line of the f -spectrum and its 
counterpart in the fiducial spectrum. 

The fiducial spectrum can be phased (effectively adjusting the acquisition phase) so that its constituent lines all 
appear upright. Thus, the answer to the N = 3 Deutsch problem may be determined by inspecting the /-spectrum 
for inverted lines. Each balanced function produces at least one inversion. For constant functions all lines are upright. 
This provides a solution state NMR scheme for conclusively answering the N = 3 Deutsch problem with just one 
application of the evolution stage to the thermal equilibrium input state. 



IV. PULSE SEQUENCE CONSTRUCTION: EXPERIMENT 

A. Molecular parameters and coupling architecture 

A saturated solution of 13 C labeled alanine in D2O provided the qubits (see Fig. ||). We denote the carboxyl carbon, 
qubit 2, the a carbon, 1 and the methyl carbon, 0. Protons were decoupled using a standard heteronuclear decoupling 
technique. Scalar couplings are J21 = 56 Hz, J10 = 36 Hz and J20 = 1.57 Hz. For J20, the linewidth is comparable 
to the relevant interpeak splitting and thus the coupling strength differs from the interpeak distance. However, the 
SWAP method of gate construction allows for a creation of an antiphase doublet between these peaks without any 
prior knowledge of J20, giving an accurate method for calculating its value (see appendix A). The T\ relaxation times 
were determined via an inversion relaxation sequence and T2's by a single spin echo. Results are 7i(2) = 20.3s, 
Ti(l) = 2.82s, Ti(0) = 1.45s, T 2 (2) = 1.25s, T 2 (l) = 0.417s, and T 2 (0) = 0.702s where the argument labels the 
qubit. 



qubit 

H 3 C^ 



NH 2 

FIG. 6. Qubit representation in the alanine molecule. 

Quadratic term gates between the pair of qubits 1 and 2 or the pair and 1 can be implemented using the directly 
coupled realization, Eq. (|27j). Here the refocusing scheme of Eq. (|3|) suggests 

[eC ad ] = [l/4J 2 i] tot - [180°]° - [l/4J 21 ] tot - [180°]° - [90°]1, - {90°t z (37) 

and 

Kl AD ] = [l/4J 10 ] tot - [180°]* " [V4Jio] tot - [ISO ], 2 , ~ [90°]°_ 2 - {90°t z (38) 

where n represents any axis in the xy plane. Note that in each of these sequences the refocusing scheme eliminates 
the coupling between qubits and 2. Unfortunately, for C/q° ad , which is required for /g and /10, the directly coupled 
realization is inadequate since 1/2 J20 = 0.32 s which is comparable to the smallest X2. However, one of the indirectly 
coupled realizations illustrated in Fig. ||or Fig. El where i = 2, k = 1 and j = 0, provides a construction which requires 







qubit 1 c q ubit2 



CH' 



OH 
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additional pulses but is sufficiently faster. The approach based on swapping to a controlled-NOT gate (see Fig. ^) is 
preferable despite the fact that it appears to require more pulses than that based on swapping to a scalar coupling 
gate. The reason is that the initial rotation of the algorithm (see Fig. |lj) may be chosen with n = — y. Then the first 
two rotations on spins 1 and cancel, leaving the scalar coupling term between and 1 as the initial operation on 
these qubits (see Fig. ||). However, these are initially in thermal equilibrium and this operation has no effect and may 
be omitted. The resulting pulse sequence has duration 0.056 ms, which is significantly faster than the direct version 
using [l/2J2o] 20 , and comfortably within the shortest relaxation time. The pulse sequence compiles to 



[t^A.,] = [90°]° - [90^ - [l/2Jio] 10 - [90% [1/2J 21 



[90°£ - [1/2J 10 ] 10 - [90% - [90°]° - [1/2J 10 ] 10 - 

[90°]^- [90°]% - [90^ - [90°]L Z (39) 

where 

[1/2J 21 ] 21 = [l/4J 21 ] tot - [180°]° - [l/4J 21 ] tot - [180°]° (40) 

and 

[1/2J 10 ] 10 = [l/4J 10 ] tot - [180°]* - [l/4J 10 ] tot - [180<t • ( 41 ) 

Again the refocusing schemes eliminate the coupling between qubits and 2. Thus, only the linear coupling chain 
— 1 — 2 is used to carry out the algorithm's evolution stage; for all admissible functions the value of J20 is irrelevant 
during this step. 



The pulse sequences for /g are obtained by concatenating those of Eqs. (37) - (p 



B. Pulse parameterization 

The experiments were performed at room temperature using a Bruker 500-DRX spectrometer and an inverse 
detection probe. The bulk of the experimental effort was devoted to constructing [C^q°j AD ] according to Eq. ( |39| ) and 
consequently to implementing the algorithm for /g and fio- 

All realizations required qubit selective 90° rotations and most also used selective 180° rotations for refocusing. 
Shaped Gaussian pulses with 10 % truncation and durations of 0.7ms for qubits and 1 and duration 0.5ms for qubit 
2 gave sufficient selectivity. No hard pulses were used. The power parameters for selective rotations were found 
by applying the selective pulse to the thermal equilibrium state followed by signal acquisition; this procedure was 
repeated, starting at a low power and increasing this with each application. The resulting acquired signals with 
maximal and minimal intensities correspond to 90° and 180° rotations. The effectiveness of the resulting 90° rotation 
can be examined partially by applying the following sequence to the thermal equilibrium state: 

[90% - delay T - [90T co8 e - +sln e - - acquire (42) 

where T is small enough to neglect any scalar coupling effects. For a rotation through exactly 90° the acquired 
signal's intensity is proportional to cos (27wT) where v is the offset relative to the receiver frequency. Furthermore, if 
the axis of the second rotation compared to that of the receiver remains constant as T varies, then the phase of the 
acquired spectrum should remain constant. We found that this was resulting signal dependence on T, particularly in 
the region where the signal intensity approaches a minimum, is fairly sensitive to variations in the pulse power. We 
used this to perform fine adjustments to the power levels for 90° pulses. 

Construction of the real pulse sequence corresponding to Eq. (p9) is accomplished by supplying the correct parame- 
ters for the active section of the pulse sequence (such as pulse phases and duration of delays) followed by computation 
of the correct passive section parameters. The division of the pulse sequence into active and passive parts is reflected 
in the different techniques required for realizing the correct operations. 

C. Passive section parameterization 

The passive section of the algorithm's pulse sequence (see Eq. (p0|)), is implemented by adjusting qubit-dependent 
phase shifts after signal acquisition. For the fiducial case and any / the passive section parameter contains absolute 
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phase information (unlike the relative rotation axes of the active section) given by the algorithm and assorted cor- 
rections which arise from Zeeman evolution during the active section. However, the fact that the algorithm answer 
criterion is based on comparing /-spectra line phases to those in the fiducial spectrum means that for any / only 
(j)^ — 0g d are important where are the passive section parameters for / and </>g d those for the fiducial case. Such 
relative values are determined by the algorithm and are independent of the acquisition phase and choice of the arbi- 
trary initial rotation phases (provided that these remain constant for all admissible functions and the fiducial case). 
It follows that an acceptable procedure is to begin by acquiring the fiducial spectrum, phase this to give upright 
absorption-mode multiplets and save the resulting post-acquisition phase shifts. Then, for any /, the active section 
of the algorithm pulse sequence is applied followed by post-acquisition phase shifts, <j?f — 0g d , added to the values 
stored after processing of the fiducial spectrum. The resulting output spectra will be phased in accordance with those 
of Fig. |. 



D. Active section parameterization 

The active section of the algorithm (see Eq. (|l9|)) is implemented by applying a sequence of spin selective Gaussian 
shaped pulses with the appropriate rotation axes and refocused delays of the appropriate durations. These parameters 
are derived from the corresponding algorithm sequence; for example, those for fg are derived from Eq. (|3^). Applying 
a pulse program with rotation axes and delays taken literally from Eq. ( |39| ) is unsatisfactory because this would ignore 
effects such as the ongoing Zeeman evolution discussed earlier. In fact, the bulk of the effort is devoted to tuning 
these parameters correctly. 

Axes of rotation are determined by pulse phases, which can be specified with a resolution of 360°/65536 = 0.0054932° 
on the Bruker equipment. In general it is not straightforward to associate spectrometer phases with absolute rotation 



axes such as those required by Eq. (|39|). However, Eq. (24) shows that it is only the relative axes and hence relative 



pulse phases which are important. The primary difficulty in adjusting these relative phases correctly arises from 
evolution under the Zeeman terms during pulse applications and delays. There are three cases: 

i) Zeeman evolution during delays, 

ii) Zeeman evolution experienced by spectator qubits during application of any selective pulse and 

iii) Zeeman evolution experienced by the rotated qubit during a selective pulse. 

During a delay of duration t the Zeeman terms produce qubit-dependent rotations about z through angles 2i{v 3 t 
where v 3 is the offset frequency of qubit j. This is corrected by re-expressing the idealized delay term as 

[t] tot = [t] tot _ [27r ^ 2 _ [2?w4] 1 _ _ r _ 2m/lt] 2 _ 1 _ (43) 

and noting that if the Zeeman terms are correctly included then the spectrometer actually implements [i] tot — [27ri/ 2 t] ^~ 

[27r^ 1 t]^ — [27w°i]^. The corrections [— 2-Kv 2 t\ 2 z — \—2ttv 1 £^ — [— 2irv t\° z are then rearranged according to the 
compilation scheme with the result that for qubit j the axes of rotation of succeeding rotations must be shifted 
through 2-kvH and the final phase shift must be altered according to (ft 3 — + cj> 7 ~ 2-kvH. 

During a selective pulse of duration t to any qubit the same rules apply to spectator qubits. The remaining 
complication is the effect of Zeeman terms on the qubit to which the selective pulse is applied. These can be 
calculated analytically for rectangularly shaped pulses by passing to a rotating frame; the results are similar to those 
for delays |32|] . However, to the best of our knowledge, no such comparable result exists for Gaussian shaped pulses. 
Our approach is to adjust the phases experimentally by considering the form of the spectrum acquired after each 
spin selective rotation. This sometimes requires readout pulses, judiciously selected to reveal key features of the 
density operator without going to the extreme of complete state tomography. This technique depends on a complete 
knowledge of the density operator prior to acquisition. Although this is not satisfactory for general purpose quantum 
computation, it is sufficient to demonstrate the essential dynamics. Furthermore it should be possible to simulate the 
exact effects of simultaneously Zeeman evolution and rotation under Gaussian shaped pulses. 

In principle, evolution under the scalar coupling terms also provides erroneous effects during the pulse sequence. 
Fortunately, the [l/4Jy] iJ delays are exempt from such errors since these are designed to utilize these couplings. It 
remains to assess the effects during 5 /zs delays between selective pulses (required by the spectrometer for switching 
purposes) and during the selective pulses. The extent of such evolution can be estimated crudely by comparing t to 
1/2 Jy, which represents the "maximum" evolution under the relevant scalar coupling. Thus a crude measure of the 
fractional error is i/jT - = ^2Jy where t is the duration of the delay or pulse. For a 5 /is delay this yield at most 
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0.00056 and for the 0.7 ms pulses and the strongest coupling, 0.078. Thus to a fair approximation the scalar coupling 
terms can be neglected during selective pulses. 

The remaining experimental issue is to construct the refocused delays of Eqs. ( fie] ) and (^lj). The durations were 
determined experimentally by comparing the signal acquired after the delay to that calculated theoretically. For 
[1/4 Jio] ° the optimal delay was found to be 6.5 ms and for [l/4J 2 i] tot it was 3.95 ms. The respective theoretical 
values were 6.9 ms and 4.5 ms. The discrepancies can only be partially explained by coupling evolution during the 
1.4 ms required for the refocusing pulses. In principle, the phases of the two 180° pulses within any refocused 
delay relative to those elsewhere in the pulse sequence are unimportant; only the relative phases of the two matter. 
These were again determined experimentally by the general process described for the selective pulses. An additional 
improvement is attained by adding the results of two experiments in which the sequences are identical except that in 
the second the axes of the refocusing pulses were both offset by 180°. This is done using the spectrometer's phase 
cycling routines. In theory the pulse sequences are identical; in practice the refocusing appears to be more accurate 
using this scheme. 



E. Spectral output and results 
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FIG. 7. 13 C output spectra for alanine: (a) Fiducial spectrum, (b) fi spectrum, (c)/42 spectrum, (d) /3 spectrum, (e) f$ 
spectrum, (f) /s spectrum, (g) /e spectrum, (h) /V spectrum, (i) fg spectrum, (j) fg spectrum, and (k) /io spectrum. In each 
spectrum the leftmost multiplet corresponds to qubit 2, the central to qubit 1 and the rightmost to qubit 0. Insets provide 
enlargements of all antiphase and doubly antiphase multiplets. 



For each representative function listed in Table | signal acquisition takes place immediately after implementation of 
Iff. Figure (^provides selected experimental spectra that are phased so that —I 3 X product operator terms correspond to 
upright multiplets. In general it is necessary to adjust the post-acquisition phase for each multiplet. The Bruker-500 
spectrometer allows only two adjustments, a constant and linear phase shift, per processed spectrum. The entire 
spectrum can be displayed with the correct phases by inserting, immediately before acquisition, a qubit selective 180° 
rotation about the appropriate axis in the xy plane for the remaining multiplet. If the axis of rotation is chosen 
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correctly relative to the density operator terms for the qubit then it has the same effect as a qubit-dependent phase 
shift. This requires knowledge of the final density operator and is not useful for general purpose quantum computation. 
In our experiment it is merely used for convenient display of the /-spectra. The arrangements of in-phase, antiphase 
and doubly antiphase multiplets agree with those predicted from pj, indicating that the correct dynamical evolution 
for the algorithm had been implemented. An interesting feature is that the doubly antiphase multiplets for qubits 
and 2 are well resolved in the fg and /io cases (see Fig. 0(e)). For qubit 2 the small J20 as well as the larger 
J21 splitting are clearly discernable. Similarly for qubit 0. The improvement in resolution of the J20 splitting in 
comparison to the inphase case stems from the fact that the resolution of an antiphase doublet is always better than 
that of an inphase doublet whenever the linewidth is comparable to the splitting (see appendix A). It follows that 
this indirectly coupled gate realization offers a method for determining the values of small couplings. 

An estimate of errors for the most complicated case, fg, may be obtained by applying a selective 90° readout pulse 
about the x axis immediately after Uf. Ideally the readout spin multiplet should remain while the others disappear. 
The average amplitudes of the residual signals for the latter lie between 14% and 31% of the average amplitude of 
the corresponding lines with no readout. The ability to extract the Deutsch problem solution for N < 3 via pure 
phase information in the output spectrum, in contrast to amplitude information, has clearly mitigated such errors. 
The most likely source of error are imperfections in the selective rotations. The most complicated case, fg, required 
eight selective 90° rotations, during each of which possible effects of scalar coupling evolution during application were 
ignored. Indeed, for the indirectly coupled realization of U 20 the total duration of all eight selective rotations, 5.6 ms, 
is 63% of 1 /2 J 2 i = 8.9 ms. To the best of our knowledge there is no satisfactory analytical understanding of the effects 
of coupling dynamics during selective pulses. At the very least such effects need to be modeled in terms of typical 
quantum information processing operations and must be applicable to shaped selective pulses. 

This work was supported, in part, by the DARPA and the ONR. We would also like to thank Gary Sanders for 
useful discussion. 



V. APPENDIX A 

Lineshapes are predicted via x and y magnetization components, which are solutions of the Bloch equations, ]37| 

M X (A.) = ^M,—^-^— (44) 



1 

1 + TfT$u%~+ T|(Aw) 2 



My(Auj) = u> x T 2 M - ■ ^ , , ^ A , A2 (45) 



where Aw is the offset from the resonance frequency, uji = 7B1 where B\ is the perturbing field strength, Mq the 
magnetization after relaxation and T\ and T2 are the relaxation times. For weak perturbations T 2 T 2 lu 2 <C 1 



2A 

l + 4(Auj/Auj h y< 

AAuj/ajh 
1 + 4(Auj/Auj h y 



where A — wxT^Mq is the maximum amplitude attained by either magnetization and Au>h — 2/T2 is the half width 
of My, By setting 

u := Aw/Aw A (48) 

and 

m x , y := M x<y /A (49) 

then 

m *W = ( 5 °) 

m > } = rrb (51) 

describe the lineshapes in frequency units of u>h and magnetization units of A. Suppose that two lines with resonance 
(angular) frequencies separated by 2irJ are superimposed. In units of linewidth the separation is 
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s := (2irJ)/Auj h . (52) 
Thus an inphase doublet in absorption mode is described by 

2u 2u 

m[u) = l + 4( U - S /2) 2 + l + 4( U + s/2) 2 (53) 

It is straightforward to show that this has three distinct extrema whenever s > l/y/3 and only one when s < 



In the former case two maxima (peaks) of height l/2s 2 (^/l + 1/s 2 — 1) are separated by y 2s^/l + s 2 — (1 + s 2 ) and 
are equidistant from a minimum (interpeak trough) of height 2/(1 + s 2 ). In the latter case the two formerly distinct 
peaks collapse into a single peak of height 2/(1 + s 2 ). 

On the other hand an antiphase doublet in absorption mode is described by 

2u 2u 
m[u) = l + 4( M -s/2) 2 ~ l + 4(u + s/2) 2 (M) 

For all s > this has two extrema; a maximum and minimum of equal amplitude and which are separated by 



2\Js 2 - 1 + v 7 -? 4 + s 2 + 1/(2^3). (55) 

Thus for < s < l/y/3 the two peaks of the antiphase doublets will, in principle, always be distinct whereas those of 
the inphase doublet always collapse into a single peak. 

It may seem that this would provide a method for determining J based on the peak separation. However, unless 
Au>h is known, this is impossible as Eq. ( |55| ) only applies when the separation is given in units of linewidth and the 
spectrometer only supplies such measurements in units of Hz. Nor would it be possible to extract such information 
from the peak heights as these require A, another typically unknown variable. 

A convenient way around these difficulties is to phase the antiphase doublet so that the lines are in dispersion mode. 
Here the spectrum is described by 

/ \ _ 2u~s 2u + s 

m[U) ~ l + 4(u-s/2) 2 l + 4( U + S /2) 2 (&bj 

The number and arrangement of extrema depends on s but for s < \/3, which is the region of most interest, there are 
two maxima of equal intensity equally spaced from a single minimum. The maxima and minima are 

(57) 



1 + s 2 

The two roots of m are situated between the maxima and are separated by 



2(i + Vi+l 5 ) 

23 j- (58) 



Auo = Vl + s 2 . (59) 
The two maxima, which are easily identifiable, are separated by 



Atw - V s2 + 1 + 2 ^1 + s 2 - (60) 

The ratio 

™ 1 i „2 

(61) 



4(i + VTT7 2 ) 

is independent of the units in which the peak amplitudes are measured and can be determined from the spectrometer 
output. Inversion of Eq. ( pl| ) equation yields s. Then, Eq. ( |60| ) gives Au max (in units of linewidth). This is related to 
the value determined from the output spectrum (in units of Hz) by a factor of Au>h- This gives the linewidth, Auj^, 
which together with s gives J. It is important to note that all that this method requires is the ability to create an 
antiphase doublet. For weak scalar couplings it is possible to do this without any knowledge of J provided that a 
network of strong couplings connects the relevant spins. 
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